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ABSTRACT 

Surface  reconstruction  from  unorganized  sample  points  is  an 
important  problem  in  computer  graphics,  computer  aided 
design,  medical  imaging  and  solid  modeling.  Recently  few 
algorithms  have  been  developed  that  have  theoretical  guar¬ 
antee  of  computing  a  topologically  correct  and  geometrically 
close  surface  under  certain  condition  on  sampling  density. 
Unfortunately,  this  sampling  condition  is  not  always  met  in 
practice  due  to  noise,  non-smoothness  or  simply  due  to  in¬ 
adequate  sampling.  This  leads  to  undesired  holes  and  other 
artifacts  in  the  output  surface.  Certain  CAD  applications 
such  as  creating  a  prototype  from  a  model  boundary  require 
a  water-tight  surface,  i.e.,  no  hole  should  be  allowed  in  the 
surface.  In  this  paper  we  describe  a  simple  algorithm  called 
Tight  Cocone  that  works  on  an  initial  mesh  generated  by 
a  popular  surface  reconstruction  algorithm  and  fills  up  all 
holes  to  reconstruct  a  water-tight  surface.  In  doing  so,  it 
does  not  introduce  any  extra  point  and  produces  a  trian¬ 
gulated  surface  interpolating  the  input  sample  points.  In 
support  of  our  method  we  present  experimental  results  with 
a  number  of  difficult  data  sets. 
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1.  INTRODUCTION 

The  problem  of  approximating  a  surface  in  three  dimen¬ 
sions  from  its  point  samples  appears  in  many  applications  of 
science  and  engineering.  Various  formulations  of  the  prob¬ 
lem  are  possible  with  different  requirements  on  the  input 
and  output.  Here  we  focus  on  the  problem  of  computing  a 
piecewise  linear  surface  that  approximates  the  original  sur¬ 
face  from  which  only  a  set  of  discrete  sample  points  are  given 
as  input. 

Because  of  its  widespread  application,  this  problem  has 
been  studied  intensely  in  recent  years.  A  very  early  paper  on 
the  problem  was  by  Boissonnat  [6]  who  proposed  a  sculpting 
of  the  Delaunay  triangulation  for  reconstruction.  A  more  re¬ 
fined  sculpting  strategy  was  designed  by  Edelsbrunner  and 
Mucke  [14]  in  their  a-shape  algorithm.  Bajaj,  Bernardini 
and  Xu  [4]  used  a-shapes  for  reconstructing  scalar  fields  and 
three  dimensional  CAD  models.  The  a-shapes  require  a  uni¬ 
form  sampling  and  an  appropriate  choice  of  a  according  to 
the  sampling  density.  Hoppe  et.  al  [18]  reconstructed  the 
surface  by  using  the  zero  level  set  of  a  distance  function 
defined  by  the  input  point  sample.  Curless  and  Levoy  [11] 
extended  this  distance  function  approach  to  multiple  range 
scans.  The  reconstruction  quality  was  further  improved  by 
Turk  and  Levoy  [20].  This  approach  produces  nice  results 
though  the  reconstructed  surface  may  not  interpolate  the 
input  data  points  and  it  requires  the  normal  information 
at  the  sample  points  in  the  input.  Bernardini  et  al.  [5] 
designed  the  ball-pivoting  algorithm  that  also  requires  the 
normals  in  the  input.  Gopi,  Krishnan  and  Silva  [17]  pro¬ 
posed  to  project  the  sample  points  with  their  neighbors  on 
a  plane  and  then  lift  the  local  two  dimensional  Delaunay 
triangulations  to  reconstruct  the  surface.  Projections  were 
also  used  by  Kobbelt  and  Botsch  [19]  where  they  took  the 
advantage  of  the  hardware. 

The  first  algorithm  that  has  been  proposed  with  theoreti¬ 
cal  guarantees  in  reconstruction  is  due  to  Amenta,  Bern  and 
Kamvysselis  [1].  This  algorithm  called  Crust  exploits  the 
structures  of  the  Voronoi  diagram  of  the  input  point  set  to 
reconstruct  the  surface.  Amenta,  Choi,  Dey  and  Leekha  [2] 
introduced  the  Cocone  algorithm  which  improved  Crust 
both  in  theory  and  practice.  Boissonnat  and  Cazals  [7]  de¬ 
signed  another  algorithm  based  on  natural  neighbors  fol¬ 
lowing  the  development  of  Crust.  Funke  and  Ramos  im¬ 
proved  the  theoretical  complexity  of  the  Cocone  algorithm 
[15].  Recently  Cohen-Steiner  and  Da  have  designed  another 
Delaunay  based  surface  reconstruction  algorithm  [10].  All 
these  algorithms  work  nicely  on  dense  data  sets,  but  they 
face  difficulty  if  the  data  contains  undersampling.  In  most 


cases  these  algorithms  produces  holes  or  other  artifacts  in 
the  vicinity  of  the  undersampling.  Many  applications  such 
as  CAD  need  water-tight  surfaces  so  that  the  surface  bounds 
a  solid.  Amenta,  Choi  and  Kolluri  [3]  devised  Power  Crust 
to  meet  this  goal.  However,  this  algorithm  introduces  many 
extra  points  in  the  output  and  also  does  not  produce  a  tri¬ 
angulated  surface.  Edelsbrunner  [13]  proposed  a  Morse  the¬ 
oretic  approach  to  reconstruct  surfaces  which  also  produces 
water-tight  surfaces.  However,  the  algorithm  may  not  re¬ 
cover  some  small  features  of  the  surface  even  if  it  is  densely 
sampled.  Along  the  same  line  Giesen  and  John  [16]  designed 
another  algorithm  based  on  the  flow  complex.  This  algo¬ 
rithm  may  not  interpolate  the  sample  points  and  introduces 
extra  points  in  the  output. 

The  purpose  of  this  paper  is  to  announce  a  simple  algo¬ 
rithm  called  Tight  Cocone  that  can  reconstruct  surfaces 
that  are  water-tight.  The  algorithm  first  computes  a  prelim¬ 
inary  surface  using  Cocone  which  almost  completes  the  re¬ 
construction  except  at  the  vicinity  of  the  undersamplings.  A 
subsequent  marking  and  peeling  phase  completes  the  recon¬ 
struction  by  filling  all  holes.  This  phase  computes  the  out¬ 
put  surface  as  the  boundary  of  a  union  of  Delaunay  tetrahe- 
dra.  This  guarantees  that  the  surface  cannot  have  any  hole 
though  it  may  contain  nonmanifold  properties  where  under¬ 
sampling  is  extreme.  In  achieving  water-tight  ness  Tight 
Cocone  does  not  introduce  any  extra  point  and  finds  tri¬ 
angles  from  the  Delaunay  triangulation  to  fill  up  the  gaps 
produced  by  Cocone. 

2.  PRELIMINARY  RECONSTRUCTION 

The  input  to  our  algorithm  is  a  point  set  in  three  dimen¬ 
sions.  We  assume  that  this  point  set  has  been  sampled  from 
a  compact  surface  5d3  that  has  no  boundary.  We  use  a 
Voronoi  based  approach  to  obtain  an  initial  approximation 
to  S  from  its  point  sample.  Amenta  and  Bern  [1]  showed 
that  the  Voronoi  cells  are  long  and  thin  along  the  direction 
of  the  normals  at  each  sample  point  if  the  sample  is  suffi¬ 
ciently  dense.  Based  on  this  important  observation  they  pi¬ 
oneered  the  surface  reconstruction  algorithm  called  Crust. 
In  a  subsequent  work,  Amenta,  Choi,  Dey  and  Leekha  pro¬ 
posed  the  Cocone  algorithm  that  improved  Crust  both  in 
theory  and  practice.  Further  improvements  were  carried  out 
by  Dey  and  Giesen  [12]  which  include  a  detection  of  under¬ 
sampling  in  the  input  data.  We  implemented  this  modified 
Cocone  algorithm  which  we  use  to  obtain  a  preliminary 
surface  out  of  the  input  sample.  We  review  this  algorithm 
briefly  below;  a  complete  description  can  be  found  in  [12]. 

2.1  Definitions 

Voronoi  diagram  and  its  dual  Delaunay  triangulation  con¬ 
stitute  the  main  data  structure  in  our  algorithm.  Let  P  be 
a  finite  set  of  points  in  M3.  The  Voronoi  cell  of  p  €  P  is 
given  as  Vp  =  {x  G  M3  :  \/q  E  P  —  { p },  \\x  —  p\\  <  ||x  —  q\\}. 
The  sets  Vp  are  convex  polyhedra.  Closed  faces  shared  by 
two  and  three  Voronoi  cells  are  called  Voronoi  facets  and 
Voronoi  edges  respectively.  The  points  shared  by  four  or 
more  Voronoi  cells  are  called  Voronoi  vertices.  The  Voronoi 
diagram  Vp  of  P  is  the  collection  of  all  Voronoi  cells,  faces, 
edges  and  vertices.  It  defines  a  cell  decomposition  of  M3 . 

The  Delaunay  triangulation  Dp  of  a  set  of  points  P  is 
dual  to  the  Voronoi  diagram  of  P.  The  convex  hull  of  four 
or  more  points  in  P  defines  a  Delaunay  cell  if  the  intersec¬ 
tion  of  the  corresponding  Voronoi  cells  is  not  empty  and 


there  exists  no  superset  of  points  in  P  with  the  same  prop¬ 
erty.  Analogously,  the  convex  hull  of  k  <  3  points  defines  a 
( k  —  1) -dimensional  Delaunay  face  if  the  intersection  of  their 
corresponding  Voronoi  cells  is  not  empty.  Every  point  in  P 
is  called  a  Delaunay  vertex.  The  collection  of  Delaunay  cells 
and  their  faces  defines  a  decomposition  of  the  convex  hull  of 
all  points  in  P.  This  decomposition  is  a  triangulation  where 
the  Delaunay  cells  are  tetrahedra  if  the  points  are  in  general 
position. 

In  our  case  P  is  the  input  point  sample  of  a  surface  S  C 
M3.  It  turns  out  that  the  Voronoi  cells  are  elongated  along 
the  normals  to  the  surface  S  if  P  is  sufficiently  dense  for 
S.  This  fact  is  used  by  Amenta,  Bern  and  Kamvysselis  to 
estimate  the  normals  at  the  sample  points  with  poles  [1]. 

Poles:  The  farthest  Voronoi  vertex  p+  in  Vp  is  called  the 
positive  pole  of  p.  The  negative  pole  of  p  is  the  farthest 
point  p~  G  Vp  from  p  such  that  the  two  vectors  from  p  to  p+ 
and  p~  make  an  angle  more  than  j .  We  call  vp  =  p+  —  p, 
the  pole  vector  for  p.  If  Vp  is  unbounded,  p+  is  taken  at 
infinity,  and  the  direction  of  vp  is  taken  as  the  average  of 
all  directions  given  by  unbounded  Voronoi  edges. 

The  pole  vector  vp  estimates  np,  the  normal  to  S  at  p. 
Therefore,  the  plane  passing  through  p  with  the  normal  as 
vp  approximates  the  tangent  plane  at  p.  A  conservative 
estimate  of  this  plane  is  given  by  thickening  it  around  p. 
This  motivates  the  following  definition  of  cocone  which  turns 
out  to  be  useful  for  the  surface  reconstruction  [2,  12]. 

Cocone :  The  set  Cp  =  {y  £  Vp  :  Z((y  —  p),  vp)  >  is 
called  the  cocone  of  p  where  Z((p—  p),  vp)  denotes  the  acute 
angle  between  the  supporting  lines  of  the  vectors  y  —  p  and 
Vp.  In  words,  Cp  is  the  complement  of  a  double  cone  (clipped 
within  Vp)  centered  at  p  with  an  opening  angle  ^  around 
the  axis  aligned  with  vp.  See  Figure  1  for  an  example  of  a 
cocone. 


Figure  1:  A  Voronoi  cell  Vp  is  elongated  along  the 
normal  np.  The  pole  vector  p+  —  p  approximates  np. 
The  cocone  Cp  is  the  region  in  Vp  between  the  two 
cones  at  p  (right). 


2.2  Cocone  algorithm 

The  Cocone  algorithm  proceeds  as  follows.  Each  sample 
chooses  a  set  of  triangles  from  the  Delaunay  triangulation  of 
the  sample  P  whose  dual  Voronoi  edges  are  intersected  by 
the  cocones  defined  at  the  sample.  All  such  chosen  triangles 
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Figure  2:  Preliminary  surface  with  COCONE  may  have  holes.  The  triangles  bordering  the  holes  are  shaded 
darker. 


over  all  samples  are  called  the  candidate  triangles.  If  the 
sampling  density  is  sufficiently  high,  these  candidate  trian¬ 
gles  lie  close  to  the  original  surface  S.  Further  they  have 
normals  oriented  nearly  in  the  same  direction  as  those  at  its 
three  sample  vertices  [2].  A  subsequent  manifold  extraction 
step  extracts  a  manifold  surface  out  of  this  set  of  candidate 
triangles.  This  manifold  is  homeomorphic  and  geometrically 
close  to  S.  All  these  claims  are  theoretically  proved  [2]. 

This  algorithm  works  nicely  when  the  input  point  set  P 
samples  S  densely.  However,  in  practice  the  input  often  un¬ 
dersamples  the  surface.  This  undersampling  may  be  caused 
by  non-smoothness,  inadequate  sampling  or  noise.  It  is  our 
experience  that  the  algorithm  computes  many  undesirable 
triangles  near  undersampled  regions.  The  appearance  of  the 
undesirable  triangles  can  be  attributed  to  the  fact  that  the 
samples  in  the  undersampled  regions  do  not  have  the  nor¬ 
mals  reliably  estimated  and  as  a  result  they  choose  ‘garbage’ 
triangles  as  candidate. 

In  a  subsequent  work  Dey  and  Giesen  [12]  proposed  a 
method  to  detect  the  undersampled  regions  from  the  input 
point  set.  The  algorithm  which  is  called  Boundary  detects 
the  sample  points  that  lie  in  the  undersampled  regions.  The 
Voronoi  cells  of  these  sample  points  are  not  long  and  thin 
along  the  normals  to  the  surface.  This  is  detected  by  two 
conditions  called  ratio  and  normal  conditions  respectively. 
The  ratio  condition  tests  the  ‘skinniness’  of  the  Voronoi  cells 
while  the  normal  condition  tests  if  its  elongation  matches 
with  those  of  its  cocone  neighbors.  We  refer  to  [12]  for  de¬ 
tails. 


The  Cocone  algorithm  is  modified  to  take  advantage  of 
the  detection  of  sample  points  in  the  undersampled  regions. 
It  calls  Boundary  and  then  lets  only  those  sample  points 
to  choose  their  candidate  triangles  which  are  not  marked 
by  Boundary.  This  modification  removes  most  of  the  un¬ 
wanted  triangles  near  the  undersampled  regions  creating 
holes  in  the  surface.  In  Figure  2  we  show  reconstructions 
with  this  Cocone  algorithm. 

3.  WATER-TIGHT  RECONSTRUCTION 

The  modified  Cocone  algorithm  detects  undersampled 
regions  and  leaves  holes  in  the  surface  near  the  vicinity  of 
undersampling.  Although  this  may  be  desirable  for  recon¬ 
structing  surfaces  with  boundaries,  many  applications  such 
as  CAD  designs  require  that  the  output  surface  be  water¬ 
tight ,  i.e.  a  surface  that  bounds  a  solid.  Formally,  we  define 
water-tight  surface  as  a  2-complex  embedded  in  M3  whose 
underlying  space  is  same  as  the  boundary  of  the  closure  of 
a  3-manifold  in  M3.  Towards  this  goal  we  design  a  very 
simple  but  elegant  algorithm  Tight  Cocone  to  compute 
water-tight  surfaces  from  an  input  point  sample. 

The  overall  idea  of  Tight  Cocone  is  to  designate  the 
Delaunay  tetrahedra  computed  from  the  input  sample  as  in 
or  out  according  to  an  initial  approximation  of  the  surface 
and  then  peeling  off  all  out  tetrahedra.  This  leaves  the  in 
tetrahedra,  the  boundary  of  whose  union  is  output  as  the 
water-tight  surface.  The  output  of  Cocone  is  taken  as  the 
initial  approximated  surface  possibly  with  holes  and  other 
artifacts. 


By  being  the  boundary  of  the  union  of  a  set  of  tetrahedra, 
the  output  surface  has  to  be  water-tight.  However,  it  is  not 
only  the  water-tightness,  but  also  the  geometric  proximity 
to  the  original  surface  that  are  desired  for  the  output.  For 
this  we  make  some  decisions  in  the  algorithm  based  on  the 
following  principle. 

Principle  of  locality.  The  undersampling  is  mostly  local. 

This  means  that  Cocone  computes  most  of  the  intended 
surface  except  with  holes  that  are  locally  repairable.  Of 
course,  if  this  principle  is  not  obeyed,  we  will  still  get  a 
water-tight  surface  but  it  may  not  be  close  to  the  original 
one. 

3.1  Marking 

The  Cocone  algorithm  as  described  in  the  previous  sec¬ 
tion  computes  a  preliminary  surface  possibly  with  holes  and 
other  artifacts  at  the  undersampled  regions.  The  sample 
points  in  the  rest  of  the  surface  have  their  neighborhoods 
well  approximated.  Specifically,  the  set  of  triangles  incident 
to  these  points  form  a  topological  disk  We  call  the  points 
whose  incident  triangles  form  a  topological  disk  good.  The 
rest  of  the  points  whose  incident  triangles  do  not  form  a 
topological  disc  are  called  poor. 

Definition  1.  The  union  of  triangles  incident  to  a  good 
point  p  is  called  its  Umbrella  denoted  as  Up. 

The  marking  of  tetrahedra  walks  through  the  Delaunay 
triangulation  in  a  depth  first  manner  using  the  vertex  and 
triangle  adjacencies.  It  maintains  a  stack  of  pairs  (p,  a) 
where  p  is  a  good  point  and  a  is  a  tetrahedron  incident  to 
p  marked  out.  Suppose  the  pair  (p,  a)  is  currently  popped 
out  from  the  stack.  The  umbreall  Up  locally  separates  the 
tetrahedra  incident  to  p  into  two  clusters,  one  on  each  side; 
see  Figure  3.  The  cluster  that  contains  a  is  marked  out  since 
a  is  already  marked  out.  The  other  cluster  gets  the  mark¬ 
ing  in.  This  is  done  by  initiating  a  local  walk  from  a  that 
traverses  all  tetrahedra  through  triangle  adjacency  without 
ever  crossing  a  triangle  in  Up  and  marking  each  tetrahedron 
as  out.  The  rest  of  the  tetrahedra  that  are  not  encountered 
in  this  walk  get  the  in  marking.  During  this  local  walk  in 
the  out  cluster,  when  a  vertex  q  of  Up  is  reached  through  a 
tetrahedron  cr\  the  pair  (q,cr')  is  put  into  the  stack  if  q  is 
good  and  is  not  explored  yet,  see  Figure  3. 

Now  we  face  the  question  of  initiating  the  stack.  For  this 
we  assume  that  Dp  is  augmented  with  ‘infinite’  tetrahedra 
that  are  incident  to  a  triangle  on  the  convex  hull  of  P  and  a 
point  at  infinity.  The  stack  is  initiated  with  a  good  point  on 
the  convex  hull  paired  with  an  incident  infinite  tetrahedron. 

For  most  of  the  data  in  practice  the  surface  computed 
by  Cocone  is  well  connected,  i.e.  all  triangles  incident  to 
good  points  can  be  reached  from  any  other  good  point  via 
a  series  of  triangle  adjacencies.  Assuming  this  connectivity 
of  the  preliminary  surface  computed  by  Cocone,  the  above 
procedure  marks  all  tetrahedra  that  are  incident  to  at  least 
one  good  sample  point.  However,  the  tetrahedra  all  of  whose 
vertices  are  poor  points  are  not  marked  by  this  step.  We  call 
them  poor  tetrahedra. 

The  poor  tetrahedra  whose  vertices  lie  in  a  single  under¬ 
sampled  region  tend  to  be  small  due  to  the  principle  of  lo¬ 
cality.  We  can  choose  to  mark  them  either  way  and  the 
peeling  process  later  takes  out  or  keeps  them  in  according 
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Figure  3:  The  umbrella  of  p  has  three  triangles 
pqr:prs  and  pqs.  This  umbrella  separates  the  the 
tetrahdera  incident  to  p  into  two  clusters,  the  up¬ 
per  cluster  {abspi  asqp ,  abpq ,  bqrp ,  brsp}  and  the  lower 
cluster  {cqrp:csrp,cqsp}.  Suppose  the  walk  entered  p 
with  the  pair  ( p,bprs ).  While  marking  all  tetrahedra 
in  the  upper  cluster  as  out ,  the  pair  (q:  bprq)  is  en¬ 
tered  into  the  stack  since  q  is  a  good  and  unexplored 
point  (right). 

to  their  marking.  In  any  case,  the  surface  gets  repaired  in 
the  undersampled  region.  For  example,  see  Figure  4. 


Figure  4:  The  four  vertices  marked  with  dark  cir¬ 
cles  border  a  hole.  The  poor  tetrahedron  with  this 
four  vertices  is  marked  in.  The  sharp  tip  p  of  the 
long  tetrahedron  is  a  good  point  which  marks  it  out. 
When  this  tetrahedron  is  peeled,  the  triangle  oppo¬ 
site  to  p  fills  the  hole  partially.  The  other  triangle 
of  the  hole  also  gets  into  the  output  by  a  similar 
peeling. 

Other  poor  tetrahedra  that  connect  vertices  from  differ¬ 
ent  undersampled  regions  tend  to  be  big.  If  such  a  big  poor 
tetrahedron  lies  outside  the  intended  surface,  we  need  to 
take  it  out.  So,  it  should  be  marked  out.  On  the  other  hand 
if  this  big  poor  tetrahedron  lies  inside  the  intended  surface, 
we  need  to  mark  it  as  in.  Otherwise,  a  large  void/tunnel  in 
the  surface  is  created  by  taking  out  this  tetrahedron.  We 
eliminate  this  dilemma  using  the  principle  of  locality.  If 
a  poor  tetrahedron  has  a  triangle  with  vertices  from  the 
same  undersampled  region,  then  that  triangle  is  small.  The 
poor  tetrahedra  lying  inside  the  intended  surface  have  to 
be  reached  by  the  peeling  process  that  peels  away  all  out 
marked  tetrahedra.  This  means  that  the  inner  poor  tetra¬ 
hedra  have  to  be  reached  through  a  small  triangle.  We  take 
this  observation  into  account  during  peeling  while  dealing 
with  the  poor  tetrahedra  and  defer  designating  them  during 


Figure  5:  The  boundary  of  the  union  of  peeled  tetrahedra  as  peeling  process  progresses. 


the  marking  step. 

3.2  Peeling 

After  the  marking  step,  a  walk  is  initiated  to  peel  off  tetra¬ 
hedra  that  are  marked  out  and  some  others.  The  boundary 
of  the  union  of  the  remaining  tetrahedra  form  the  water 
tight  surface.  This  is  also  the  surface  of  the  union  of  the 
peeled  tetrahedra,  see  Figure  5. 

The  walk  maintains  a  stack  of  surface  triangles  that  form 
the  boundary  of  the  union  of  the  tetrahedra  peeled  so  far. 
It  is  initiated  with  all  convex  hull  triangles.  At  any  generic 
step,  a  triangle,  say  £,  is  popped  out  from  the  stack.  One  of 
the  tetrahedra  incident  to  t  is  already  peeled.  If  the  other 
incident  tetrahedron,  say  a  is  marked  in ,  we  put  t  in  the 
output  list.  Otherwise,  there  are  two  possibilities  depending 
on  whether  a  has  been  already  peeled  or  not.  If  a  is  already 
peeled,  the  triangle  t  separates  two  out  tetrahedra  and  thus 
should  not  be  part  of  the  output.  In  case  a  has  not  been 
peeled,  the  walk  should  move  past  t  into  a  if  a  is  not  a  poor 
tetrahedron.  In  this  case  we  replace  t  with  the  other  three 
triangles  of  a  into  the  stack.  In  case  a  is  a  poor  tetrahedron, 
the  walk  is  allowed  to  move  past  t  only  if  t  is  not  the  smallest 
triangle  in  a.  This  is  done  to  protect  peeling  of  the  inner 
poor  tetrahedra  as  we  discussed  before.  Notice  that  if  a  is 
a  poor  tetrahedron  outside  the  intended  surface,  it  will  be 
eventually  reached  by  the  peeling  process  at  triangles  other 
than  the  smallest  one.  But,  if  a  is  a  poor  tetrahedron  inside, 
it  can  only  be  reached  from  outside  through  its  smallest 
triangle  due  to  the  principle  of  locality. 

The  walk  terminates  with  the  surface  triangles  in  the  out¬ 
put  list  when  there  is  no  more  triangle  to  process  from  the 
stack. 

4.  EXPERIMENTAL  RESULTS 
4.1  Examples 

In  Figure  7  we  show  the  results  of  our  implementation  of 
Tight  Cocone  on  some  difficult  data  sets.  In  Mannequin 
there  are  undersamplings  in  eyes,  lips  and  ears  which  pro¬ 
duce  holes.  Tight  Cocone  closes  all  these  holes.  In  partic¬ 
ular,  in  the  ear  there  is  a  relatively  large  hole  since  points 
cannot  be  sampled  for  occlusion.  This  hole  is  nicely  filled. 
The  Pig  data  has  severe  undersampling  in  the  hoofs,  ears 
and  nose.  They  are  mostly  due  to  the  fact  that  these  thin 
and  highly  curved  regions  should  have  more  sample  points 
to  capture  the  features  properly.  Tight  Cocone  fills  all 
holes  and  produces  a  water-tight  surface  for  this  difficult 
data  set.  The  machine  part  Connector  has  undersam¬ 
pling  mainly  due  to  non-smoothness.  Sharp  edges  and  cor¬ 
ners  cause  undersampling  which  cannot  be  avoided  by  any 


finite  sampling.  All  anomalies  caused  by  this  non-smooth 
regions  are  repaired  to  produce  a  water-tight  surface.  The 
undersampling  in  the  highly  curved  regions  such  as  the  neck 
in  Dinosaur  is  also  gracefully  handled  by  Tight  Cocone. 

The  code  for  Tight  Cocone  has  been  released  for  public 
use.  It  can  be  downloaded  from  [9]. 

4.2  Timings 

The  time  and  space  complexities  of  Tight  Cocone  are 
dominated  by  those  of  the  three  dimensional  Voronoi  dia¬ 
gram  computation.  Thus,  for  an  input  sample  of  n  points, 
the  algorithm  runs  in  0(n2)  time  and  space  in  the  worst 
case.  However,  in  practice  we  do  not  observe  this  quadratic 
behavior.  In  Table  6  we  show  the  breakup  of  the  timings  of 
Tight  Cocone  for  a  number  of  data  sets  on  a  PC  with  733 
Mhz  Pentium  III  CPU  and  512  MB  memory.  The  code  was 
compiled  with  g++  compiler  and  at  the  01  level  of  optimiza¬ 
tion.  CGAL  2.3  is  used  for  all  computations.  The  timing 
is  broken  into  the  three  major  steps  of  Tight  Cocone,  the 
Delaunay  triangulation,  the  Cocone  surface  generation  and 
the  final  postprocessing  to  make  the  surface  water-tight.  It 
is  clear  from  the  table  that  the  Delaunay  triangulation  is  the 
most  time  consuming  step.  All  computations  are  done  with 
filtered  floating  point  arithmetic  for  numerical  robustness.  If 
we  are  willing  to  compromise  the  robustness,  floating  point 
arithmetic  can  be  used  which  usually  provides  a  gain  of  a 
factor  of  two  in  computing  time. 


object 

#  points 

Delaunay 

time(sec.) 

Cocone 

time(sec.) 

Mark/peel 

time(sec.) 

Cactus 

3337 

6.56 

2.72 

1.38 

Pig 

3511 

2.82 

2.7 

1.47 

Mechpart 

4102 

3.68 

3.27 

1.73 

Earcanal 

8459 

16.54 

7.16 

3.22 

Cat 

10000 

8.51 

8.22 

4 

Knot 

10000 

10.01 

11.18 

5.29 

Mannequin 

12772 

9.31 

10.08 

4.89 

Dinosaur 

14050 

12.63 

12.16 

5.91 

Fandisk 

16475 

16.16 

12.44 

6.3 

Club 

16864 

23.12 

13.47 

6.34 

Connector 

26793 

41.58 

18.5 

9.68 

Female 

30432 

38.56 

27.61 

13.28 

OlLPUMP 

30936 

28.68 

23.92 

11.88 

Heart 

37912 

39.99 

30.36 

14.5 

Horse 

48485 

90.24 

42.07 

19.45 

Brevi 

56152 

54.39 

41.82 

20.28 

Hand 

74315 

123 

64.21 

29.83 

Figure  6:  Time  data. 


Water-tight  Mannequin 


Water-tight  Pig 


Holes  in  high  curvature  regions 


Water-tight  Connector 


Anomalies  near  sharp  edges  and  corners 


Water-tight  Dinosaur  Holes  in  high  curvature  regions 


Holes  are  filled 


Figure  7:  Results  of  Tight  Cocone.  Second  column  shows  the  holes  in  the  preliminary  surface  while  the  third 
column  shows  how  they  are  filled. 


5.  CONCLUSION 

We  designed  a  simple  and  elegant  algorithm  to  produce 
water-tight  surfaces  out  of  point  samples.  The  algorithm 
guarantees  water- tight  ness.  It  does  not  introduce  any  extra 
point  to  do  so.  The  output  is  always  a  subcomplex  of  the 
Delaunay  triangulation.  Thus,  one  can  use  the  Delaunay 
tetrahedra  inside  the  water-tight  surface  to  produce  a  three 
dimensional  triangulation  of  the  model.  Our  experiments 
show  the  effectiveness  and  the  efficiency  of  the  algorithm  on 
many  example  data  sets. 

In  some  applications  the  output  surface  needs  to  produce 
certain  intended  boundaries.  The  Cocone  algorithm  is  de¬ 
signed  for  that  purpose.  Although  it  detects  the  intended 
boundaries,  it  also  produces  holes  where  undersampling  hap¬ 
pens  unintentionally.  Is  it  possible  to  recognize  the  intended 
boundaries  while  filling  up  other  holes?  The  problem  is 
ill-posed,  but  is  there  a  proper  formulation  of  the  problem 
which  may  lead  to  an  algorithm  with  certain  reasonable  as¬ 
sumptions  on  input?  We  plan  to  pursue  this  line  of  research. 
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